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ABSTRACT 

We discuss theoretical predictions concerning the evolution of globular cluster Pop. II stars 
vis-a-vis current estimates of standard errors in the determination of nuclear burning rates. 
Numerical evaluations are given for the dependence of TO and HB luminosity from the rate 
of the relevant nuclear reactions. We conclude that evolutionary predictions appear rather 
solid in this respect, with a maximum 3a error of about 1 Gyr in the evaluation of cluster ages 
derived from the calibration of the difference in luminosity between HB and TO. However, 
current evaluation of the original He content, as given on the basis of the R parameter, 
should wait for a much better determination of the C12(a, 7)016 reaction before reaching a 
satisfactory accuracy. 



1 INTRODUCTION 

The theoretical evolutionary scenario for Pop. II globu- 
lar cluster stars has long been recognized as a key tool to 
reach a close insight into the early evolution of the Galaxy 
and, in turn, on the history of the Universe. In recent times, 
theoretical constraints concerning the age of galactic globu- 
lar clusters have been widely debated in the literature, be- 
cause of a possible contradiction between those ages and 
the age of the Universe as derived from current estimates 
of the Hubble constant (see, e.g., van den Bergh, 1994). In 
this context, of particular relevance is a recent paper by 
Chaboyer (1995) who discussed the degree of confidence to 
be assigned to theory, investigating globular cluster ages for 
suitable variations in the adopted input physics. 

In this paper we will follow a similar approach, dis- 
cussing the evolutionary behaviour of low mass Pop. II stars 
on the basis of current uncertainties in the nuclear burning 
rates relevant for similar structures. The consequences of an 
overall variation of H burning rates on the cluster Turn Off 
(TO) luminosity has been already discussed in the above 
paper by Chaboyer. Here we will enter in more details, 
showing that the contribution of pp or CNO burning rates 
has to be separately taken into account, since an increase 
(or a decrease) of both pp and CNO rates has opposite and 
compensating effects on TO luminosity. Moreover, we will 
discuss the influence of the relevant reaction rates on current 
evaluation of Horizontal Branch (HB) structures, investigat- 
ing in particular HB luminosities and lifetimes. 

By inserting in such a scenario suitable values of stan- 
dard errors for nuclear cross sections we will finally evaluate 
the corresponding errors for selected evolutionary features, 
namely the TO luminosity, the difference in luminosity be- 
tween HB and TO and the parameter R constraining the 
amount (Yrj) of original He in cluster stars. We will find 



that current uncertainties in the nuclear burning rates have 
only a minor influence on the theoretical predictions con- 
straining the cluster age, but a not negligible influence on 
current estimates of Yo- 

2 H BURNING PHASES 

All the evolutionary computations presented in this paper 
have been produced by adopting the FRANEC stellar evo- 
lutionary code, as already described in previous papers (see, 
e.g., Chieffi & Straniero, 1989 and reference therein). 

To investigate H burning evolutionary phases a suitable 
set of stellar models with masses in the range 0.7 to 0.9 Mq 
has been followed with the standard version of the code from 
the Zero Age Main Sequence (ZAMS) phase till the onset of 
the He flash. We assumed Y=0.23 and Z=0.0001 or 0.001 
as suitable original chemical composition covering the range 
of metallicities found in globular cluster stars. The compu- 
tations have been thus repeated for different assumptions 
about the astrophysical factors governing the efficiency of 
the relevant burning rates. 

Figure 1 compares the canonical evolutionary track cov- 
ering the H burning phases of a 0.8 Mq (Z=0.0001) model 
with similar tracks but with the pp cross section ( 

(fpp) in- 
creased by a factor 1.3 (left panel) or the CNO cross section 
(acNo) increased by a factor of 3 (right panel). 

A closer insight into the role of nuclear burning rates in 
the evolution is given in Table 1 where we report selected 
structural parameters of the 0.8 Mq model for the labeled 
choices about the efficiency of reactions. Left to right one 
finds the luminosity of the track TO (logL^o), the evolu- 
tionary time at the TO (tTo)i the luminosity of the Red 
Giant (RG) clump (logL c j ump ), the luminosity at the onset 
of the He flash (logLy; as ^), the mass of the He core at this 
time (Mc) and the Red Giant Branch (RGB) evolutionary 
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Figure 1. The evolutionary path in the HR diagram of standard stellar models (full lines) compared with similar models but with pp 
burning rates increased by a factor 1.3 (left panel) or with CNO rates increased by a factor of 3 (right panel). 



Table 1. Selected evolutionary parameters for a 0.8 Mq star 
with Z=0.0001, Y=0.23 and for the labeled assumptions on nu- 
clear reaction rates 

Rates logL TO t TO logL c/Mmp logLj. lash Mc t rgb 

Gyr Myr 

Standard 0.401 13.49 2.239 3.272 0.509 80.4 
ppxl.3 0.414 13.71 

CNOx3 0.365 13.28 2.193 3.306 0.503 78.1 

time (t RGB-defined as the time the star elapses between 
log(L/L ) =1.7 and the He flash). 

Data in table 1 show that increased pp or CNO rates 
have opposite effects on the luminosity of the evolutionary 
models and, in particular, on the luminosity of the track 
TO. Such an occurrence can be easily understood on a quite 
simple basis: let us here only recall that the track TO is 
marked by the onset of CNO burning, and this onset is ob- 
viously trought forward when CNO rates are increased or 
pp rates are decreased. 

As expected, one finds that variations in pp rates do not 
affect at all the evolution of RG structures. On the contrary, 
increasing the CNO rates moderately decreases the luminos- 
ity of the RG bump marking the encounter of the H burning 
shell with the chemical discontinuity produced by the first 
dredge up. Moreover, one finds that variations modify the 
RG core-luminosity relation, since the He flash is reached 
at a larger luminosity (AlogL fi ash =0.034) but with a lower 
mass of the He core (AMc=-0.006). However, evolutionary 



times along the RG branch are practically unaffected. 

The set of evolutionary tracks has been finally used to 
produce cluster isochrones under the various assumptions 
about nuclear reaction rates. One finds that the sensitiv- 
ity of cluster TO luminosity to variations in burning rates 
is marginally dependent on the mass of the evolving stars 
and thus on the assumed cluster age. As a whole, one can 
safely assume AlogLxo~ -0.03 as the maximum variation 
expected for a decrease of pp rates by 30% or, alternatively, 
for an increase of CNO rates of 200% . In terms of incre- 
ments one thus finds; 

A logL TO - O.lSpp - 0.0158CNO 

where Spp and 5CNO represent the fractional variation 
in the nuclear burning rates. 

3 HE BURNING PHASE 

Since the luminosity of the isochrone TO is the clock mark- 
ing the cluster ages, the results discussed in the previous 
section will be easily translated into theoretical uncertain- 
ties of age estimates. However, experimental (i.e. observa- 
tional) constraints on such a luminosity can be derived from 
the observed TO magnitudes only if preliminary estimates 
of the cluster distance modulus (DM) are available and suit- 
able evaluations for the bolometric corrections are provided. 
The DM is often -but not in all cases- evaluated using HB 
stars as standard (theoretical) candles; in this way one can 
directly calibrate the difference in luminosity between HB 
and TO (logL hb-To) m terms of cluster ages (see, e.g., 
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Figure 2. Standard ZAHB sequences (dashed-dotted lines) for Z=0.0001 or Z=0.001 compared with ZAHB sequences computed by 
(only) increasing by a factor 3 the CNO reaction rates (short dashed lines), and by also taking into account the predicted decrease of 
Me (full lines). Dotted lines connect selected models with the same mass. 



Renzini & Fusi Pecci 1988). According to such a procedure, 
the HB luminosity behaves as a further ingredient entering 
in the age determination, whose dependence on the assumed 
burning rates will be investigated in this section. 

HB luminosity plays a role also in the estimates of an- 
other relevant evolutionary parameter, the amount of orig- 
inal He (Yo). Following a well established procedure sug- 
gested by Iben (1968), the estimate of Yo relies on the ob- 
servational value of the parameter R, the ratio between the 
number of HB stars to the number of RG stars above the 
HB luminosity level (see, e.g., Renzini & Fusi Pecci 1988). 
Theoretical calibration of this parameter, besides depend- 
ing on HB luminosity, relies then on evolutionary predic- 
tions concerning stellar lifetimes both in the RG and in the 
HB phases. Evolutionary times along the RGB have been 
already discussed. In the following we will thus focus our 
attention also on HB lifetimes, to reach an indication of the 
accuracy of theoretical estimates for both cluster ages and 
original He content. 

3.1 H burning reactions 

As already discussed, pp reaction rates have no influence 
on the mass of the He core at the onset of the He flash 
(Mc). Nor do variations in these rates affect the distribu- 
tion of chemical elements throughout the structure of a new- 
born HB star. Bearing in mind that H shell burning in HB 
stars is dominated by CNO reactions, one easily foresees 



that HB evolution should be largely unaffected by changes 
in pp rates. This prediction is indeed confirmed by numeri- 
cal experiments. By decreasing the pp rate by 30% one finds 
that the Zero Age Horizontal Branch (ZAHB) locus remains 
unchanged, with only a slight increase of the masses located 
at a given effective temperature. As an example, with the 
quoted decreased rate, the mass of the model at logTe=3.85 
is M=0.802 Mq against M=0.796 M© for the canonical case. 
At the same time He burning evolution of the models does 
not show relevant variations. 

CNO burning rates deserve more attention, since they 
affect both Mc and the energy sources in HB stars. The 
rather complex behaviour of ZAHB models with CNO rates 
increased by a factor of 3 is shown in Figure 2, where canoni- 
cal ZAHBs for Z=0.0001 or Z=0.001 are compared with sim- 
ilar loci but computed i) by increasing only CNO rates and, 
ii) by considering also the predicted decrease of Mc. The 
increase of the CNO rates drives a corresponding increase of 
the ZAHB luminosity, whereas the decrease of Mc works in 
the opposite direction. However, Figure 2 shows the unex- 
pected result that the net affect of these two contributions 
to ZAHB luminosity substantially depends on the assumed 
metallicity. In more detail, at logTe=3.85, one finds that the 
luminosity of the ZAHB decreases by about AlogL=0.017 
when Z=0.0001, whereas it increases by about the same 
amount if Z=0.001. The same Figure 2 shows that such a 
behavior follows the larger variation in luminosity induced 
by the increased CNO rates in the Z=0.001 case, due to 
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Figure 3. Dependence of HB lifetimes on the C12(a, 7)016 reaction rate 



the greater efficiency of CNO burning in similar moderately 
metal rich structures. According to such evidence, one finds 
that when Z=0.0001 the decrease of Mc succeeds in pushing 
the final ZAHB luminosity below the original value, whereas 
for Z= 0.001 the increase in luminosity prevails over the ef- 
fect of Mc. 

3.2 He burning reaction 

The 3a reaction starts affecting HB structures from the final 
phase of RG evolution, governing the onset of the He flash 
and, in turn, the mass of the He core in a new born ZAHB 
star. Numerical experiments show that decreasing the 3a 
rate by a factor 0.7 increases Mc by AMc= 0.003. According 
to the pioneering paper by Sweigart & Gross (1976) we know 
that HB luminosity depends on Mc according to the relation 
AlogL£f5=3.4A Mc. Consequently the quoted variation of 
Mc should, alone, induce an increase of 0.010 in the ZAHB 
luminosity. Numerical experiments disclose that such a pre- 
diction is only slightly increased by the contribution given 
by decreased 3a rates. As a result, one finds that for the 
given variation of 3a rates ZAHB luminosity will increase 
by about AlogL=0.012. In the meantime, the lifetime tjjb 
for the entire phase of central He burning increases by about 
0.7% . 

To study the role of the C12(a,7)016 reaction we com- 
puted selected HB models by increasing that rate (CTC12) 
by selected factors. As expected, HB lifetime increases 
since the completion of the chain of reactions is favored 



Table 2. The expected variation of the labeled evolutionary 
parameters when the various nuclear cross section are increased 
by a factor of 2. 



App ACNO 
AlogL TO 0.1 -0.015 
AlogL HB ±0.009 



A3a AC* 

-0.040 
-0.02 



12 



0.10 



and more energy is obtained from the He burning. Fig- 
ure 3 shows that HB lifetimes do not depend linearly on 
°"C12- This could be expected since HB lifetimes cannot be 
affected by (Ten m the limits where all or no carbon is con- 
verted in oxygen. Figure 3 however shows that in the range 
°"C12/ ,<J ci2 — 2 °ne can adopt a linear approximation given 
by At hb /t hb ~ O.lOAi7ci2/0'ci2 ■ Note tnat one nnds 
HB lifetimes increased by values larger than predicted by 
previous estimates reported in Renzini & Fusi Pecci (1988). 

Table 2 summarizes the results of all the numerical ex- 
periments, giving the dependence of TO luminosity, HB lu- 
minosity and HB lifetimes on changes in the rates of the 
relevant nuclear reactions. 



4 STANDARD ERRORS 

We are now in the position to link numerical experiments 
discussed in previous sections with current estimates for 
standard errors in burning rates in order to estimates stan- 
dard errors in the predicted evolutionary features. 
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Table 3. The 3<r errors for the labeled reactions and the cor- 
responding uncertainties in logL^Q, logL^fg and logL^fg_-pQ. 
Last row gives the total 3<x error obtained by quadratically com- 
pounding the pp, CNO and 3a uncertainties 



Rates Errors AlogL^Q AlogL^fjj AlogL^f g_ 



pp 5% 

CNO 51% 

3a 45% 
3o-tot 



0.005 
0.008 

0.009 



0.005 
0.018 
0.019 



Z=10~ 4 
0.005 
0.003 
0.018 
0.019 



Z=10~ 3 
0.005 
0.013 
0.018 
0.023 



For pp reactions, from the recent European coopera- 
tion on nuclear data (NACRE97) one derives S(0) = 3.94 (1 
-0 025)10~ 25 Mev barns (DegPInnocenti, private communi- 
cation). As a consequence we assume in the following a stan- 
dard error of 5% (at 3 a). Concerning CNO rates, Rolf & 
Rodney (1988) give for the the main reaction N14(p,7)015 
a range of value as large as S(0)= 2-10 kev barns. How- 
ever, this range appears greatly reduced in Schroder et al. 
(1987), which gives a 1 a error of about 17%. In the fol- 
lowing we will adopt this last result. According to Rolf & 
Rodney (1988) the error in the 3a reaction can be again esti- 
mated to be of the order of 15%. However the uncertainty in 
the C12(tv, 7)016 reaction rate is much worse. According to 
Caughlan & Fowler (1988) reaction rates usually adopted in 
current evolutionary computations could be wrong by a fac- 
tor of 2 up or down. In a more recent determination Buch- 
mann (1996) limits S(300) to the range 62-270 keV barns. 

By relying on the above estimates it is now possible to 
carry out the exercise of deriving expected standard errors 
on the luminosities of TO and HB. The results of this ex- 
ercise are reported in Table 3, where we list for each given 
reaction the assumed 3a errors and the corresponding ex- 
pected variation in evolutionary parameters. 

Concerning the calibration of logL^o in terms of clus- 
ter age, since dlogL^o~dlogtg (see, e.g., Castellani & 
Degl'Innocenti 1995), for cluster age of 14 billion years we 
obtain a standard error of the order of 0.3 Gyr, in reason- 
able agreement with the values given by Chaboyer (1995). 
From the same table one finds that the HB luminosity is 
substantially affected only by CNO and 3a rates. 

In the last two columns of table 3 we report the uncer- 
tainties of the difference between the TO and the HB lumi- 
nosities. In the evaluation of this quantity it is necessary, 
in principle, to consider separately the two cases Z=0.0001 
and Z=0.001, due to the opposite behaviour of the HB lumi- 
nosities discussed in section 3.1. However the general result 
is that the uncertainty in logL hb~TO is> a t most, of order 
0.023 .As a results, when calibrating ages through the dif- 
ference in magnitude between HB and TO, the error is of 
the order of magnitude of 0.8 Gyr. 

Data in Table 2 allow us finally to investigate the influ- 
ence of " nuclear" errors on current estimates of the amount 
of the original He, as based on the calibration of the R pa- 
rameter in terms of Yo . By recalling that R is theoretically 
defined as the ratio of He burning lifetime to RG lifetime 
at luminosity larger than the HB luminosity, one finds that 
the relevant reactions are CNO and 3a which affect the HB 
luminosity only, and C12+a which governs HB lifetime. 

By relying on available estimates of the lifetime along 
the RG branch (see, e.g., Bono et al. 1995) one finds that 



the uncertainty in the HB luminosity level translates into an 
uncertainty of 3.5% in the evolutionary time of an RG above 
this level, and a corresponding error in the parameter R. On 
the basis of theoretical calibration of such a parameter (see 
again Bono et al. 1995) one concludes that the contribution 
to the Yo error is rather negligible, less than 0.01. 

Since current calibrations of R are based on Caughlan 
et al. (1985) reaction rates an error of the order of the fac- 
tor 2 has been assumed to investigate the relevance of the 
C12(a, 7)016 reaction rates uncertainties. With such an er- 
ror the C12(a, 7)016 reaction affects the He burning evolu- 
tionary time by about 10%. Again on the basis of current 
R calibration, it turns out that the theoretical calibration of 
He is affected by an error of order AY$ ~ 0.02 which appears 
now of some relevance. As an example, increasing &C12 by 
a factor two would move Yq from Yq=0.25 to Yq=0.23. 



5 CONCLUSIONS 

In this paper we have discussed the reliability of the- 
oretical prescriptions concerning the evolution of globular 
cluster Pop. II stars vis-a-vis current estimates for standard 
errors in the determination of nuclear burning rates. We 
computed stellar evolutionary models adopting different as- 
sumptions for the reaction rates involved during the H and 
He burning phases, evaluating the dependence of selected 
evolutionary features on the nuclear reaction rates. The 
main results can be summarized as follow: 

i) The uncertainties in pp rates sligthly affect only the TO 
luminosity; 

ii) The CNO rates uncertainties may introduce a small vari- 
ation in the TO and HB luminosities. Due to the different 
efficiency of CNO burning, the effect of a given variation of 
CNO on HB luminosity depends on the assumed metallicity; 

iii) An increase (or a decrease) in the pp rate has the oppo- 
site effect on TO luminosities of an increase (or a decrease) 
in the CNO rates, almost removing any sizeable variation in 
logL TO ; 

iv) The uncertainty of the 3a rates affects the HB luminosity 
but has a negligible effect on the HB evolutionary time; 

v) The uncertainty of the C12(a,7)016 rate seems to intro- 
duce a non-negligible uncertainty in the HB lifetime. 

On this basis we investigated the relevance of the cur- 
rent nuclear rection rates uncertainties to theoretical cali- 
brations concerning cluster ages and the original He content 
Yo- We conclude that cluster ages derived from the calibra- 
tion of the difference in luminosities between HB and TO are 
affected by maximum 3<r error of about 1 Gyr. However, it 
appears that the theoretical calibration of the parameter R 
should wait for better determination of the C12(a,7)016 
reaction before reaching a satisfactory accuracy. 

We would like to thank S. Degl'Innocenti for useful dis- 
cussions and comments. 
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